Time evolution of the chiral phase transition during a spherical expansion 
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We examine the non-equilibrium time evolution of the hadronic plasma produced in a relativistic 
\£) , heavy ion collision, assuming a spherical expansion into the vacuum. We study the 0(4) linear 

■ sigma model to leading order in a large- N expansion. Starting at a temperature above the phase 
ON I transition, the system expands and cools, finally settling into the broken symmetry vacuum state. 

We consider the proper time evolution of the effective pion mass, the order parameter (er), and the 
' particle number distribution. We examine several different initial conditions and look for instabilities 

Oh, (exponentially growing long wavelength modes) which can lead to the formation of disoriented chiral 

■ condensates (DCCs). We find that instabilities exist for proper times which are less than 3 fm/c. 
We also show that an experimental signature of domain growth is an increase in the low momentum 
spectrum of outgoing pions when compared to an expansion in thermal equilibrium. In comparison 
to particle production during a longitudinal expansion, we find that in a spherical expansion the 

^/-^ | system reaches the "out" regime much faster and more particles get produced. However the size of 

. the unstable region, which is related to the domain size of DCCs, is not enhanced. 
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O ■ I. INTRODUCTION 

m ■ 

There have been many recent investigations into the formation of disoriented chiral condensates (DCCs) following 
a relativistic heavy-ion collision . The original motivation for studying this problem was the Centauro events , 
""^5 ■ rare cosmic ray events in which a deficit of neutral pions was observed [Q . In a recent work []5| , the time evolution of 
the hadronic plasma produced in such a collision was studied by using the 0(4) linear sigma model in a longitudinal 
expansion. The large-iV expansion was used to incorporate non-equilibrium and quantum effects into the problem. 
' After performing numerical simulations to solve the time-dependent equations of motion, instabilities were found to 
exist for only a short time, and thus no significant amount of pion domains would be formed. In this work, we study 
the same problem using a spherical expansion, since at late times the expansion becomes spherical. This situation 
. , produces the most rapid cooling of the system. We would like to see if the formation of instabilities in this geometry 

■ is more pronounced than in a longitudinal expansion. 

There are two questions which should be examined. First, we want to know which types of initial conditions lead 
to the formation of instabilities in the system; and second, if instabilities do form, we want to find out if the size of 
the unstable region is large enough to see significant domain growth. To answer the first question, we examine the 
proper-time evolution of the system, starting a short time after the phase transition, where the linear sigma model 
is appropriate. We look for the effective mass of the pion to go negative during the time evolution. This signifies 
the onset of growth in long-wavelength modes, which is believed to lead to the formation of DCCs. When we have 
a case where instabilities form, we then compute the momentum distribution of outgoing pions, and compare to a 
hydrodynamical model calculation, assuming local thermal equilibrium. We find a noticeable enhancement of low 
momentum modes as compared to the hydrodynamical model. This provides an experimental signature which can be 
measured. The implication is that the system is evolving out of thermal equilibrium, which is a necessary condition 
to have significant growth of low momentum modes. 
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We prepare the initial state of the system in local thermal equilibrium, and study the evolution using scale invariant 
kinematics (v = r/t) to model the cooling of the plasma. Scale invariant kinematics are appropriate for a high energy 
spherical expansion starting from a point source MM, and imply that mean- field expectation values only depend on 



the proper time r = yi? — r 2 . We incorporate non-equilibrium and quantum effects through the use of the large- TV 
expansion. 

The paper is organiz ed a s follows. In Sec. [njwe describe the linear sigma model to leading order in a 1/N expansion. 



We then discuss in Sec. Ill our choice of coordinates and derive equations of motion for the system. We examine issues 
of renormalization and choice of suitable initial conditions. We also derive an expression for the phase space number 
density that will be used to calculate momentum distributions. In Sec. IV we discuss the calculation of the energy- 
momentum tensor, and other thermodynamic relations. In Sec. |v| we present numerical results of our simulation. In 
Appendix [A] we look at properties of the radial functions used in the expansion for the quantum modes. In Appendix 
M we discuss the transformation of the number density to physically measurable variables. 



II. LINEAR SIGMA MODEL 



The Lagrangian density for the linear sigma model in a generalized curvilinear coordinate system is given by: 

-^(,)-, 2 ] 2 }, (2.1) 



where the mesons are in a 0(4) vector, $ = (a, if). The factor of \/—g{x) : where —.9(2;) = det[g fJlU (x)}, has been 
introduced to make the Lagrangian a scalar density. The potential here is the "Mexican hat" , with degenerate minima 
at any values of $ such that $^ — v 2 . We will remove this symmetry by introducing a non-zero current term in the a 
direction. In this work, we use the convention of an implied sum over a repeated index i, which runs from 1 to N. (Here 
N = 4.) The counting for the large-A expansion is implemented by introducing a composite field \ — A(<&f — v 2 ). 
That is, we add to the Lagrangian a term given by J8|: 

[ x -\(<S> 2 (x)-v*)} 2 /4\. 

This gives an equivalent Lagrangian, 
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-^4x + ^ 2 }. (^) 

We can write the action as: 

S[$i,rJi] = J d 4 x [£[§i,x\ + ^gjx)^} . (2.3) 
We consider the generating functional, given by the path integral, 

Z\ji] = / d[ X ] / d[*,] exp[iS[$i,x;ji}] = ejq>\iW\ji]] . 

The large-iV approximation is equivalent to integrating out the <!>i variables, then performing the remaining \ 
integral using the method of steepest descent. We then Legendre transform to find the effective action, 



T[4H,X\ = W\ji] - I d 4 xy/~g(x)Mx)Mx) 

where 



2 



We then obtain, to lowest order (in large N) 



1 



(2.4) 



T[<f>i,X] = J d^x^^^g^xW^W^] 

1 .2 V 2 1 2 

iN 

+ — ^Gq (x,x;x) 
where 

(□ + X {x))G a (x, x>- x) = i5\x - x')/^gjx). 
r[0i, x) is the classical action plus the trace- log term. 

III. COOLING MECHANISM 
A. Coordinate System 

At late times following a heavy ion collision, the energy flow becomes three-dimensional. If we assume the entire 
flow is spherical, then the system can be described in terms of the fluid variables: 



r = \/ 1 2 — r 2 , 
i] = tanh _1 (r/<) = - In 



t + r 
2 ~ [t-r 

where t = t cosh 77 and r — rsinhTy. We restrict the range of these variables to the forward light cone, < t < 00, 
and < 77 < 00. The variables r and 77 are useful to describe a free spherical expansion of a plasma into the vacuum, 
with the velocity of the fluid identified as v — r/t = tanh?7 when we are at high energies so that the expansion can 
be thought of as coming from a point source (scaling limit). Minkowski's line element is given by: 

ds 2 = dr 2 - r 2 {dT/ 2 + sinh 2 77 dO 2 + sinh 2 77 sin 2 Ode/) 2 } , 

from which we can read off the metric tensor 

9nv = diag(l, — t 2 , — t 2 sinh 2 77, — t 2 sinh 2 77 sin 2 9) , 
g — t 3 sinh 2 77 sin 9 . 

This can be compared to the Robertson- Walker metric for spherical geometry, given by the line element: 

ds 2 = dr 2 

- a 2 (t) { dT/ 2 + sinh 2 77 d9 2 + sinh 2 77 sin 2 9 d(j> 2 } . 

Thus the case we consider here corresponds to a cosmological model with a fixed uniform expansion proportional to 
the proper time r, and zero curvature. 

B. Equations of motion 



We can derive the equations of motion from the effective action, (2.4). Varying the action with respect to 4>i and 
X gives: 

ip+x{x))<t>i{x)=j i (x)=H5 i0 



(□ + x(x))G (x, x') = iS 4 (x - x')/^g[x) , (3.1) 
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and the gap equation 



X (x) = \{-v 2 + 4>1{x) + NG (x,x)} 



(3.2) 



In order to give the pions mass, it is only necessary to have a current in the zero (a) direction, so that jo(x) = H 
constant. 

We determine the parameters in the model by considering the vacuum sector. Spatial homogeneity then gives: 



Xo<J = H , 

Xo = -Av 2 + Ao- 2 + XN 



k 2 dk 



2^ 2 2^Tx~ 

The PC AC (partial conservation of axial vector current) condition gives: 

d^(x) = H n {x) , 

where the axial vector current is given by 



which leads to 



A?(x) = [iri{x)d^a(x) - a{x)d f "TT l {x)} 



H = Urn 2 



Since we have defined the vacuum by xo^o = m <ro = H, we then find that <ro = f^- The coupling constant A is 
determined by the low-energy n — n scattering data, as described in B . 



We now specialize to the case when fa and x are functions of r only. We can see that (3.1) is also th e eq uation for a 
free scalar quantum field with a time-dependent mass x( r )i which is self-consistently determined by ( |3.2| ). Therefore 
we can introduce a quantum field — fa + fa . The equations for <3?i are: 



l_d_ 

t 3 dr 



d 



T"— I +X(r) }fa(r) =HS i0 



where the 4- vector x = (r, rj, 9, 



(n + x (r))fa(x) = 0, 
Then for Go we find: 

G (x,x') = (T c {fax,r), fax',r')}) 



(3.3) 



where T c corresponds to a r-ordered product following the closed-time-path formalism of Schwinger. When 
(it 1 ) — 0, this is the true Green's function. 

Following Parker and Fulling fin], we expand fa into a complete set of states, 



ds ^ [O'i.slm Ips {r)y s l m (?/, 4>) 

I in 

+ h.c] 



(3.4) 



with 



and where -0, 7r, and Y satisfy: 



y s i m (v, 0, fa) = ir s i(r))Yi m {0, < 



(3.5) 



1 d ( , d \ s 2 + 1 



3 dr \ dr 



X(r) 



4> s = o 



1 d ( . a d 
q — sinn 77 — — 



sinh r]drj\ 



dr) 



1 - 



sinh rj 



7T S 1 = 



4 



1 d 
^6de{ Sm " 89 



1 



8 2 



sin 9 V d(j)- 



+ 1(1 + 1) 



Y lm = 0. 



Here, Yi 



s 



are the usual spherical harmonics, and TT s i(r]) are a complete set of functions, discussed in Appendix 



In order to satisfy the canonical commutation relations, we require ip s (r) to satisfy the Wronskian condition: 

r s (r)Mr)-MrW s (r) = -i/r 3 , 
where the "dot" means differentiation with respect to t. Therefore we find 

[<Pi,slm{j),(j) j s , llml {T)] = i5ijS(s' - s)5 U i5 mm </T 3 , 



and 



We are now in a position to calculate (4>f). We shall choose the (Heisenberg) state vector such that the bilinear forms 
of creation and destruction operators are diagonal: 

(aj,s>i> m >al slm ) = (n s + 1) 6ij S(s' - s) 6w5 mm ' , 

(aj,s'l'm'ai,slm) — Ps 5ij S(s' — s) SwS mm ' , 
( a L'*'mA>im) = P*s 6 V S ( S ' ~ S ) S U' S rnm> ■ 

Here n s and p s are the particle and pair densities. They will be taken to be a function of s only. In addition, we will 
take n s to be a thermal distribution in the comoving frame, 

1 



s w s (r )/fc B T _ I ' 



with ut s = ^/s 2 /tq + x(to). We can choose p s — for all our simulations, since one has the freedom to make a 
Bogoliubov transformation at To so that this is true. Using the results in Appendix |a|, we then find 

ds(2n s + l)|^(r)| 2 \y si m {v,d ^)\ 2 



2tt 5 



(2n s + l)|^(r)| : 



Therefore (3.2) becomes 



X (t) =-\v 2 + \tf(T)+\N 



2tt 2 



(2n s + l)|^ s (r)| : 



(3.6) 



and is a function of r only. This completes the derivation of the equations of motion. 



C. Initial Conditions and Renormalization 



The variable r does not allow for a good WKB expansion due to the singularity at r = 0. The transformation 
U = ln(mr), where m is any mass scale (we choose m = m w ), maps the singularity to u — — oo, and allows one to 
perform a WKB expansion in the usual manner. 

Changing variables to u and rescaling the fields using the substitutions 

4> s (u) = g s (u)e- u m 3/2 
<t>i{u) = pi(u)e~ u m 3/2 , 
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we get: 



i 2 i / \ 2m / 2 

(Mr 



g a (it) = 

Pi {u) = H6 m e 3u /m 7 / 2 



(3.7) 



with 



X (u) = -\v 2 + A^m 3 e- 2 V 2 H 
' s 2 ds 



AN 



,0 2tT 2 

with s m = Ae u /m. Then g s (u) obeys the Wronskian 



(2n s + l)m 3 e^ u \g s (u)\ z , 



(3.8) 



J.W-; j 9s{u) = -lira. 



du 



du 



From (3.7), we notice that when x( u ) < 0j modes satisfying s 2 /r 2 < % grow exponentially in r. However, from ( |3.§| ), 
we see that these growing modes will quickly cause x to become positive when A is large, as is the case here. 
We can now use a WKB ansatz for g s (u): 



9s{u) 



1 



y/2mW s (u) 



exp 



where W s (u) satisfies: 



1 W s " 3 fW' 



2 W s 4\W. 



-i / W s {u')du' 
Juo 



7 TFT +^ 2 =c s 2 ( U ), 



(3.9) 



and u> s (u) = -\/s 2 + x(it)e 2u /m 2 . We will then take the initial conditions 

W s (u ) = w s (m ) 
Wi(«o)=a;i(«o), 

which correspond to the adiabatic vacuum. This allows us to introduce an interpolating number density which 
interpolates between n s (u) and the true "out" density, n ou t- 

By a WKB analysis, one can show |ll]] that Gq(x,x) has quadratic and logarithmic divergences. The quadratic 
divergence can be removed by mass renormalization. In the vacuum sector, the mass of the pion is given by Q3.g| ), 
with Xvac = ml — to 2 : 



-Aw 2 + A/ 2 + XN 



A 7.2 



k 2 dk 



o 2tt 2 2Vfc 2 + m 2 



(3.10) 



with cutoff A. We note that if we change variables in the integral to s = kr = ke u /m, (3.10) becomes: 

to 2 = -Xv 2 + A/ 2 



-\Nm 2 e~ 2u 



1 



2tt 2 2Vs 2 + e 2 « ' 

Subtracting this equation from ( |3.8| ), we obtain a logarithmically divergent expression for x : 

X ( U )=m 2 + A^(m 3 e - 2 > 2 ( M )- / 2 ) 



(3.11) 



+XN 
—XN 



s 2 ds 
2tt 2 


{2n s 


+ l)m 3 e- 2u \g s (u)\ 2 


s 2 ds 


m 2 e~ 


1 


2ir 2 


2Vs 2 + e 2u 



(3.12) 



G 



Note that the last integral is independent of u. 
The coupling constant is renormalized by taking 



1 
A 



1 



N 



k 2 dk 



X r 8tt 2 J q (fc 2 + m 2 ) 3 / 2 

i n r Sm s 2 ds 

A7 ~ 8^ J a ( a 2 + e 2 U )3/2 



One can explicitly show by using (3.13) in (3.12) that x( u ) is now completely finite. 
We will also need the value of x(u) for the initial conditions: 



x( u ) 



1 + XN 



s 2 ds 



o 2tt 
2\e- 2 ^[p i {u)p' i {u) - p1{u)] 
XNA 3 (2n Sm + 1) 
4tt 2 v /A 2 T^ 



(2n s + 1) 



4wf(u) 



+ 



A7V 



2tt 2 



(2n s + l)mV 
X(w)e 2u 



1 



2m 2 ojf(u) w s (u) 



(3.13) 



(3.14) 



D. Phase Space Interpolating Number Density 



In the expansion (3.4), we can also use a time-dependent set of creation and annhilation operators, with first order 
adiabatic mode functions, 



4>i(u 



>»?,0></>) = / ds^2[ai, s i m (u)ip°(u)y s i m (n,6, 

J lm 



h.c] 



(3.15) 



where 



and 



</>°(it) = gl{u)e-m^ 



: exp 



y/2raw s (u) 

We can define the first order adiabatic number density as: 

rig (it) = {a\{u)a s {u)) 



i J U! s (u')du' 



(3.16) 



where, for simplicity, we have suppressed the angular indices on a s and a\. 

One can show |12] that n s (u) is an adiabatic invariant, and would be the true number density in a slowly varying 
expansion. We then choose 

a s = a s (u ) 
9s(u ) = g° s (u ) 

so that the initial a and a* are the adiabatic ones. When x( u ) ~ * ™ 2 then n s (u) — ► Ti ou t, which is the true out-state 
phase space number density. 

The time-dependent creation and annihilation operators satisfy 
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<i<L n daX n . , 

-r9° s + -r9°* = o- (3.i7) 

au au 

The time-dependent operators can be connected to the time-independent operators via a Bogoliubov transformation: 

a s (u) = a(s,u)a s + (3(s,u)a_ s , (3.18) 

where a and f3 are determined by: 

a{s,u) = i[g, — g s 

\ au au 

«•■«>-'(«!&- 2?*) ■ < 3 - 19 ' 

We then find 

n s (u) = n s (u ) + \/3{s, u)\ 2 (l + 2n s (u )) . (3.20) 

Notice that at u = u$, (3(s,uq) = 0, so n s (u) = n s (uo) = n s . Since n s (uo) is the initial phase space number density, 
and at late times, becomes the out-state number density, it is an interpolating number density. 

IV. ENERGY-MOMENTUM TENSOR 

The energy-momentum tensor T^ v is defined by: 

SS = ~ [ d i x^—gT^{x)8g^{x) , (4.1) 



2 . 

with the action given by (|2.3| ) . Performing the variations, we find the "improved" energy- momentum tensor [fl3f 

+ \ [g^g^l^ - > (4-2) 



where the Lagrangian density is given by (2.2). We follow the standard practice |7|] and define the energy density and 
pressures by 

T^ v = d\&g (e, p v t 2 ,pgT 2 sinh 2 77,p0T 2 sinh 2 77sin 2 6>) . 

Next, we take expectation values of the energy-momentum tensor. One can easily show all the pressures are equal, 
P = Pv = Ps = P<p> an d that T p „ is diagonal. 

The energy-momentum tensor obeys a conservation law: 



One can easily verify that the energy conservation equation, the v = component of ( |4.3| ) takes the form: 

|i + i( e +p)=0. (4.4) 

OT T 

In thermal equilibrium, this would become the entropy conservation equation. Then we have de = Tds and e+p = Ts. 
In an ultra-relativistic fluid expansion, with Cq = |, p — |e and er 4 = constant. 
In our simulations we have verified that the local energy conservation ( |4.3| ) is valid. 
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V. NUMERICAL RESULTS 



To choose the initial conditions, we start the system at a temperature above the phase transition in thermal 
equilibrium, with all particle masses positive. The equations are solved self-consistently at the starting time to obtain 
the values of the %, (cr) and (n) fields. We fixed the value of x at the initial time as the solution of the gap equation 
in the initial thermal state. We also required that the initial expectation values of the a and if fields satisfy 

7t 2 (t ) + cr 2 (r ) = a% , 

where ot is the equilibrium value of <& at the initial temperature T. We choose T — 200 MeV, which gives <jt = 0.3 fm 
_1 . We compute the time evolution of these fields, starting at a proper time tq = 1 fm. The value of f v used in all the 
simulations is 92.5 MeV, and Xr is 7.3 (see ||). Below we show results for several sets of initial conditions. Once the 
initial values are chosen, we have the freedom to vary the first derivative of the $ field. The results with tt ^ were 
similar to those with a ^ 0, so we only show results for the latter case. We find there is a wide range of values which 
will allow the system to become unstable, 0.15 < \&\ < 4.95. This can be compared with the longitudinal expansion 
fl, where the regime of instability was much smaller, 0.25 < \&\ < 1.3. This is because the spherical expansion leads 
to a much larger negative gradient for x than the longitudinal case. 

Figures [l] and ^ show the results of the numerical simulation for the proper time evolution of the system. We 
display the auxiliary field x m units of fm~ 2 , the classical fields $j in units of fm -1 , and the proper time in units of 
fm. At proper times greater than « 10 fm, the auxiliary field reaches its vacuum value of m 2 , the a field approaches 
its vacuum value of /„., and the tt field its value of zero. This is in distinction to the longitudinal expansion, where 
even at r = 30 fm, one had not yet reached the "out" regime. 

For all of our initial conditions, the size of the unstable region (x < 0) is at most 2-3 fm. Thus we find that the 
spherical expansion does not produce larger domains than the longitudinal case. We also see that when we start the 
initial expectation value of $j in the tv\ direction that the system becomes slightly more unstable. We notice that if 
the derivative of the field is positive or zero that this is insufficient to generate instabilities. 

In Fig. |3|, we show the effect of the initial temperature on the evolution of the auxiliary field. We see that varying 
the initial temperature has little effect. In Fig. [|, we show the evolution of the x ne ld for different values of the cutoff 
A. We can see that x is independent of A, which shows that the renormalization has been carried out correctly. In 
our simulations we use the value A = 800 MeV, since A = 1 GeV is too close to the Landau pole (see ||). When one 
chooses a cutoff too close to the Landau pole the late time be havio r becomes unstable. 

Figures || and ^ show the number density calculated from (3.20), for several different proper times. Figure || is a 



case where instabilities have arisen in the system, and there is a large amount of particle production during the time 
that x nas gone negative. Figure ^| is a case with no instabilities, and there is very little particle production as a 
function of proper time. 

Figures and || show the same distributions transformed to the physical momentum p, as discussed in Appendix |^. 
The momentum p is plotted in units of m^. We compare these distributions to a hydrodynamical model calculation [see 



(B3)], where we have assumed that when the system reaches the "out" regime, the final distribution is a combination 
of a thermal distribution in the comoving frame at T c — boosted to the center of mass frame using the boost 
variable rj(r,t) (see 0). For comparison purposes, we have renormalized the thermal distributions to give the same 
center of mass energy (E = 100 GeV) as the corresponding non-thermal distributions. We see that as a result of 
the non-equilibrium evolution, there is an enhancement at low momentum independent of whether or not there are 
instabilities; however, the effect of instabilities is to greatly magnify this low momentum enhancement. 



VI. CONCLUSIONS 



In this paper we have performed numerical simulations to examine the chiral phase transition during a uniform 
spherical expansion of the hadronic plasma. We used the linear a model to leading order in a large- TV expansion and 
studied a wide range of initial conditions starting above the critical temperature for the phase transition. We found 
initial conditions which drove the system to instabilities that lead to the formation of disoriented chiral condensates. 
Due to the necessity of a rather large renormalized coupling constant, the formation of instabilities lasted only a short 



time because of the large exponentially growing quantum corrections in (3.8), and no significant amount of domain 



formation was observed. However, we find that the phase space number density for our non-equilibrium evolution 
is significantly different from one which would result from an evolution in thermal equilibrium. The experimental 
signature for domain formation is an increase in the pion particle production rate at low momentum. Our calculations 
were done in a mean-field approximation, where all the mode coupling is due to the presence of this mean field. In 
next order in large- TV, scattering in the background mean field occurs, and the possibility for re-equilibrization exists. 
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These effects will be incorporated in a future calculation. In comparison to particle production during a longitudinal 
expansion, we find that in a spherical expansion the system reaches the "out" regime much faster and more particles 
get produced. However the size of the unstable region, which is related to the size of the domain of DCCs, is not 
enhanced. 
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APPENDIX A: PROPERTIES OF THE ty S l FUNCTIONS 
In this appendix, we discuss properties of the functions n s i(r]), which are real solutions of the equation, 
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or, 




for 77 in the range: < 77 < 00. With the substitution, 



u 8 i(r))/ sinh 77 , 



we find that u s i(r)) satisfies: 




u s i = . 



The general solution is given by |T^] : 




where the normalization M s [ is given by 



M sl = V(7r/2)s 2 (s 2 + l 2 )"-(s 2 + * 2 ) • 



The completeness relation is given by: 



00 



dsTT s i{rr)TT s i(ri') = S(rj - ^/[sinhrisinhr/] . 



For the functions y s i m defined in ( |3.5| ), the addition formula is fiof] : 

X] y*lm {Vl A» 01 ) ysim 0?2, 02, 02) 




where r\ is defined by: 



cosh rj = cosh rji cosh 772 — sinh r/i sinh 772 cos 
cos 9 = cos 0i cos 02 + sin 0i s'm8 2 cos(0i — 02 ) • 
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Therefore, taking the limit, (771, 6\, </>i) — * (772, O2, (fa), or r\ — > 0, we find 



2tt 2 



By differentiating both sides of the addition formula, we can show that 



E 

I in 

E 

I in 

E 

l in 



dy slm ( v ,9,4>) 



drj 

dy s i m {v,9,(l>) 



de 



d4> 



s 

2tt 2 
2tt 2 



s 2 + 1 



s 2 + 1 



s 2 + 1 



sinh 77 . 



sinh 77 sin 2 9 . 



APPENDIX B: TRANSFORMATION TO PHYSICAL VARIABLES 

In terms of the initial distribution of particles ti s (tq) and (3 we have: 

n s (u) = n s (u ) + |/3(s,u)| 2 (l + 2n s (u )) , 

where n s (u) is the adiabatic invariant interpolating phase space number density which becomes the actual particle 
phase space number density in the comoving frame when interactions have ceased. We now need to relate this quantity 
to the physical spectra of particles measured in the lab. At late r > 77 ~ 10 fm our system relaxes to the vacuum 
and x becomes the square of the physical pion mass m 2 . The comoving center of mass energy of outgoing particles 
can then be identified with 



The actual distribution of momenta in the lab frame is a combination of the collective ( "fluid" ) motion described 
by the boost 77 from the comoving frame to center of mass frame and the comoving particle distribution. Here, the 
space-like hypersurface on which one is counting particles is at fixed proper time r/. This distribution is given by the 
Cooper-Frye formula |15| which is: 



dN dN 

E Wr E ^dv =ifMpdcT » 



(Bl) 



We identify the relativistic phase space distribution function /(cc,p) with n s (v,f). The dependence of s on the space 
time variable x and the outgoing momentum p is found from the relationship 



^( T /) = ]l^2 +7Jl2 



We choose the measured momentum p to be in the z direction 63 of our spherical coordinate system. We have that 

w'' = (cosh 77, sinh 77 e r ) 
p» = (E,pe 3 ), 

so that 

P^Ufj. = E cosh 77 — p cos 9 sinh 7/ . 
The surface on which one is counting particles is the time-like surface r = Tf with 
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Changing variables from r to ry at fixed r we then obtain 



J f(x,p)drj dcosO sinh 2 77 p p u M , (B2) 



where 

p^Ufj = E cosh rj — p cos # sinh 77 



and we have used the isotropy assumption and chosen p as the z axis. Here E = ^Jp 1 + m 2 . 

The calculation that (B2) needs to be compared with is a hydrodynamical model calculation for a local thermal 
equilibrium flow. In a hydrodynamical model of heavy ion collisions dJ^], the final spectrum of pions is given by a 
combination of the fluid flow and a local thermal equilibrium distribution in the comoving frame. One calculates this 
spectrum at the critical temperature T c (x, t) when the energy density goes below 

1 



(h/mc) 



This defines the breakup surface r c , after which the particles no longer interact so that this distibution is frozen at 
that temperature. For an ultrarelativistic gas of pions, this occurs when T c = to. The covariant form for the spectra 
of particles is again given by the Cooper- Frye formula |15|: 



E- dN 



Aitp^dp 



= J f(x,p)drj dcosd r c 3 sinh 2 77 p^u^ (B3) 



but now f(x,p) is the single particle relativistic phase space distribution function for pions in local thermal equilibrium 
at a comoving temperature T c (r c ): 

f(x,p) = {ex P KVT c ] - IK 1 . 
We have identified the left hand side of (|B3|) as n t h(p, t). 
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FIG. 1. Proper time evolution of the \ field for the following initial conditions: Solid line is for <t(to) = ot, 7r(to) = 0, and 
&{to) = — 1- Dashed line is for ff(ro) = or, 7Ti( r o) = 0, and &(to) = 1. Dotted line is for tr(ro) = or, 7Ti( r o) = 0, and &(to) = 0. 
Dot-dash line is for a(r ) = 0, 7Ti(r ) = a T , and &(r ) = -1. At T = 200 MeV, o T = 0.3 fm" 1 . 
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FIG. 2. Proper time evolution of the (a) and (7Ti) fields for the same initial conditions as Fig. [l]. 



14 




15 



1.0 




-0.5 1 1 

5 10 15 20 

x(fm) 



FIG. 4. Proper time evolution of the \ field for three different values of the cutoff A, with A = 600, 800, 1000 MeV for the 
initial conditions cr(ro) = or, vri( r o) = 0; an d <H T o) ~ — 1- 
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FIG. 5. n s (r) computed from (|3.2C|), for the initial conditions <t(to) = ctt, 7Ti(ro) = 0, and a (to) = — 1. 
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